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Abstract 

Suppose that the matrix equation AXB = C with unknown matrix 
X is given, where A, B, and C are known matrices of suitable sizes. 
The matrix nearness problem is considered over the general and least 
squares solutions of the matrix equation AXB = C when the equation 
is consistent and inconsistent, respectively. The implicit form of the best 
approximate solutions of the problems over the set of symmetric and the 
set of skew-symmetric matrices are established as well. Moreover, some 
numerical examples are given for the problems considered. 
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1 Introduction and Notations 

Let n mxn , STZ nxn : and SSTZ nxn be the set of to x n real matrices, the set of 
nxn real symmetric matrices, and the set ofnxn real skew-symmetric matrices, 
respectively. The symbols A T , A*, and \\A\\ will denote the transpose, the 
Moore-Penrose generalized inverse, and the Frobenius norm (see, for example, 
[7]), respectively, of a matrix A G lZ mxn . Further, vec (•) will stand for the vec 

operator, i.e. vec (A) = (of, aj, . . . , a^J for the matrix A = (a%, . . . , a n ) G 
K mXn , ai G R mxl , i = 1,2, . ..,n, and A ® B will stand for the Kronecker 
product of matrices A and B G lZ mxn , (see [5]). 

The well-known linear matrix equation AXB = C, where A, B, C are 
known matrices of suitable sizes and X is the matrix of unknowns, were studied 
in the case of special solution structures, e.g. symmetric, triangular or diagonal 
solution A in [TJ [21 [31 31 Q]3 [TB] using matrix decomposition such as the singular 
value decomposition (SVD), the generalized SVD, the quotient SVD, and the 
canonical correlation decomposition. In these literatures, the matrix equation 
AXB — C is consistent. But, it is rarely possible to satisfy the consistency 
condition of the matrix equation AXB = C, since the matrices A, B, and C 
occurring in practice are usually obtained from an experiment. 

An iteration method to solve the linear matrix equation AXB = C over the 
set of symmetric matrices have constructed by Peng et al. [17) . In addition, 
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Peng P15] has established an iterative method to solve the minimum Frobe- 
nius norm residual problem: min || AX'S — C\\ where X is the symmetric ma- 
trix of unknowns. Huang and Yin solved the constrained inverse eigenproblem 
and associated approximation problem for anti-Hcrmitian i?-symmetric matri- 
ces and the matrix inverse problem and its optimal approximation problem for 
i?-symmetric matrices in [5] and [5], respectively. Huang et al. gave the precise 
solutions to the minimum residual problem and the matrix nearness problem 
for symmetric matrices or skew-symmetric matrices in |10j and constructed an 
iterative method to solve the linear matrix equation AXB = C over the set of 
skew-symmetric matrices in 

This work is devoted to give the best approximate solutions of the following 
two problems, which are interesting and known as the matrix nearness problems, 
in an alternative way: 

Problem 1 For given matrices A G TZ mxn , B G W xr , and C G K mxr , let 
Sq be the set of all solutions of the consistent matrix equation 

AXB = C. 

For a given matrix Xq G TZ nxp , find X G Sa such that 



X-X Q 



min II JSC — Xq\ 

xgS g 



Problem 2 For given matrices A e TZ mxn , B G TZ pxr , and C £ TZ mxr , let 
Se be the set of all least squares solutions of the minimum residual problem 



min \\AXB - C\\ 



For a given matrix Xq G TZ nxp , find X G Se such that 



X-Xn = min \\X-X, 



xes E 



o 



In fact, the Problems [T] and [5] are to find the best approximate solution for 
a given matrix Xq G TZ nx P over the set of general solutions (Sq) and the least 
square solutions (Se) of the matrix equation AXB — C, respectively. These 
problems are known as the matrix nearness problem in the literature. The 
matrix nearness problem is very important in applied sciences and has been 
extensively studied in recent years (see, for example, [HI H21 UHl EXS])- Therefore, 
it is important to give the best approximate solutions of the problems in implicit 
forms. 

In general, numerical algorithms or iteration methods for solving these prob- 
lems are suggested in most of the works mentioned above. 

In this work, the implicit forms of the best approximate solutions to the 
problems mentioned above have been obtained over the set of symmetric and the 
set of skew-symmetric matrices using the Moore-Penrose generalized inverse. 
Moreover, some numerical examples are given via Matlab 7.5. The matrices in 
the examples have been taken from related reference work articles. 
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2 Preliminary Results 

The vector xq G lZ nxl is the best approximate solution (BAS) to the inconsistent 
system of linear equations Ax — g, where A G TZ mxn , if and only if 

1- {Ax - gf {Ax -g)> {Ax - gf {Ax - g) for all x € K nxl , 

2- x T x > x T xo for all x G lZ nxl \{x Q } satisfying {Ax — g) T {Ax — g) = 

(Axo-gf {Axo-g) 0. 

The vector xq G lZ nxl is a least squares solution (LSS) to the inconsistent 
system of linear equations Ax — g, where A G lZ nxl , if and only if 

{Ax - gf {Ax - g) > {Ax - gf {Ax - g) 

for all x G K nxl [5]. 

It is noteworthy that there may be many LSS for an inconsistent system of 
linear equations. In addition, an LSS may not be the BAS while the BAS is 
always a LSS. However, the BAS is always unique. 

We close this section by giving two auxiliary results related to the problems 
mentioned earlier and which will be used in the rest of the work. 

Lemma 1 Suppose that Sg> is the set of all solutions to the consistent system 
of linear equations Ax = g, where A G JZ mxn is a known matrix, g G "JZ mxl is 
a known vector, and x G lZ nxl is the vector of unknowns. For a given vector 
xq G lZ nxl , the vector x G S' G satisfying 

\\x - x \\ = min ||x - xo|| 

is given by 

x = A^g + {I - A^A) x - 

Proof. If x G Sc, then it can be written in the form 

x = A^g + (J - A^A) h 

for some vector h G TZ nxl [51 Theorem 6.3.2]. Thus, the problem turns into the 
problem of finding the BAS x of the system 

A*g+ (I- A*A)h = Xo (2.1) 

or equivalently the system 

(J - A 1 A) h = x Q - A^g. 
Since the matrix / — A^ A is symmetric and idempotent, it is obtained 

h=(l- A^A) x„ 

by Theorem 7.4.1 in [5]. Substituting this expression in the equation (|2.1[) . we 
get 

x = A f g + {I - A f A) x . 
So, the proof is completed. ■ 
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Lemma 2 Let Se> be the set of all least squares solutions to the system of 
linear equations Ax — g which do not need to be consistent, where A G j^/nxn 
is a known matrix, g G TZ mxl is a known vector, and x G lZ nxl is the vectors 
of unknowns. For a given vector xq G lZ n , the vector x G S' E satisfying 

\\x -x \\ = min \\x - x a \\ 

2:6 S E i 

is given by 

x = A^g + (I - A^A) xq. 

Proof. If x G Se>i then it can be written in the form 

x = A^g + (I - AU) h 

for some vector h G lZ nxl [SI Theorem 6.3.2]. This is of the same type with 
(|2.1[) . So, the remaining part of the proof can be completed easily in a similar 
way as in the proof of Lemma [TJ ■ 



It is noteworthy that the structures of x in Lemmas [T] and [2] are exactly the 
same. 



3 The Best Approximate Solutions of Problems 
1 and 2 

If it is assumed that the matrix equation AXB = C, where A € lZ mxn , B G 
K pxr , C G ji-mxr are k nown nonzero matrices and X G 7Z nxp is the matrix 
of unknowns, is inconsistent, as was the system of linear equations, then it 
may be asked to find a matrix X such that ||AX_B — C\\ is minimum, too. A 
matrix satisfying this condition is called an approximate solution to the matrix 
equation. The matrix X G 7^. nx P i s defined to be the BAS to the matrix equation 
AXB = C if and only if 



1- \\AXB-C\\ > 



AXB - C 



for all X G TZ nxp : 



2- UAH > 



X 



for all A G ft" xp \ j satisfying \\AXB - C\\ = 



AXB - C 



We note that a vector k G TZ mnxl will stand for the vector vec(K) in the 
rest of the text, where K G lZ mxn . 

It is known that the matrix equation AXB = C can be equivalently written 

as 

(B T ®A)x = c, (3.1) 

where (B T ® A) is the kronecker product (see, for detail, [50] )• Consequently, 
the solutions of a matrix equation AXB = C can be obtained considering the 
usual system of linear equations (13.11) instead of the matrix equation AXB = C. 

Now we can give the solutions of Problems [1] and [5] which are the subjects 
of the following two theorems, respectively. 
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Theorem 1 Let the matrix equation AXB = C be consistent, where A G 
TZ mxn , B G W xr , C G TZ mxr are known nonzero matrices and and X G K nxp 
is the matrix of unknowns. Then, for a given matrix Xq G lZ nxp , the matrix 
X € Sq satisfying 

Y - X 



— min \\X — Xn\ 

xes G 



is given by 



X = A^CB' +X Q - A^AXqBBK 

Proof. If X G Sg, then it can be written in the form 

X = A^CB^ + H - A^AHBB' 

for some matrix H G TZ nx P [2Jj]. Since the statement 
statement 

A 



(3.2) 

is equivalent to the 



x = (B T A)^c {B T ® A)^ (B T <g> A) 

the problem turns into the problem of finding the best approximate solution of 
the usual system of linear equations 



(B T ® A) f c (B T ® A^ (B T ® A) 



x 



or equivalently 



/ - (B T ® A)^ (B T ®A)]h = x - {B T <8 A) f c. 
Using Lemma Q] and (|3. 31) . we get 



.T 



(3.3) 



£ = (S T ®^4) t c+ [/- ( B T ® A) f (S T ® A) 
or, in the matrix form, 

X = A^CB^ + X Q - A^AX a BB\ 
So, the proof is completed. ■ 

Lemma 3 Let the matrices A, B , and C be as in Theorem [7J and assume 
that the matrix X is a least squares solution to the inconsistent matrix equation 
AXB = C . Then the matrix X can be written in the form 

X = A^CB^ + H - A^AHBB' 

for some matrix H G TZ nxp . 

Proof. The proof is immediately follows from Theorem 6.3.2 and Theorem 
7.6.3 in [5] considering the usual system of linear equations (-B T ® A)x — c 
instead of the matrix equations AXB = C because the former and the latter 
are equivalent. ■ 

Theorem 2 Let the matrix equation AXB = C be inconsistent, where A £ 
TZ mxn , B £ K pxr , C £ K mxr are known nonzero matrices and X G TZ nxp 
is the matrix of unknowns. Then, for a given matrix Xq E lZ nxp , the matrix 
X G Se satisfying 

X -X Q 



= min \\X — Xf)\ 
xeS E 



is given by 



X = A^CB^ +X - A^AXqBB^ 
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Proof. By Lemma [31 any least squares solutions of the inconsistent matrix 
equation AXB = C is in the form 

X = A^CB^ + H - A ] AHBB^ 

for some matrix H G lZ nxp . Hence, in the framework of Lemma [21 the proof is 
easily completed by proceeding as in the proof of Theorem [TJ ■ 



4 The Symmetric and Skew— Symmetric Solu- 
tions of Problems [I] and [2] 

Now, suppose that the symmetric solutions of Problems [T] and [5] are required. 
To do this, the pair of matrix equations 

AXB = C 
B T XA T = C T 

or, equivalently, the usual system of linear equations 



" B T 


<E> A ' 




Cl 


A® 


B T 


x = 


c 2 



instead of the matrix equation AXB — C is taken with c\ = vec (C) and 
C2 = vec(C T ). Then, for a given matrix Xo G SlZ nxn , in the framework 
of Theorems [1] and [5] the solution matrix X G SlZ nxn is obtained by using 



(x^j given by 



" B T 


®A ' 


t 


Cl 




' B T 


® A ' 


t 


' B T 


®A ' 


A <g 


B T 




c 2 


+ x a - 


A « 


B T 




A& 


B T 



x . 



(4.1) 



Remark 1 If the matrix Xq is not symmetric, then the matrix \ {Xq + Xq^ , 
instead of the matrix X$, is taken to find the symmetric solutions of Problems 
[7] or0 The reason for this is that the minimization problem 

minpr-X || 

is equivalent the minimization problem 



X-'-iXo 



xt 



over the subset Sg or Se in STZ nxn 
1 



since 

2 



:v--(x + x T ) 



n*-*oir 

(for example, see \H\ \T3 \j for details). 



{Xo - X T ) 



yx g sn r 



Now suppose that the skew-symmetric solutions of Problems [T] and [5] are 
required. To do this, the pair of matrix equations 



AXB = C 
B T XA T = 



-C 2 



6 



or, equivalently, the usual system of linear equations 



" B T 


<E)A ' 




C\ 


A<8 


B T 


x = 


-c 2 



instead of the matrix equation AXB = C is taken where c\ and c 2 are as 
in Remark [1] Then, for a given matrix Xq G SSlZ nxn , the solution matrix 

X 6 SSlZ nxn is obtained by using x — vec (^X^j given by 





' B T 


® A " 


t 


Cl 




' B T 


® A " 


t 




®A " 


X = 


A (g 






-c 2 


+ a;o - 


A (g 






A (g 





(4.2) 



Remark 2 If the matrix Xq is not skew-symmetric, then the matrix ^ (Xo — ) 
instead of the matrix Xq is taken to find the skew-symmetric solutions of the 
Problems[l\ or® The reason for this is that the minimization problem 

min||X-X || 

is equivalent the minimization problem 



mm 



x-\{x Q - X T ) 



over the subset Sq or Se in SSlZ nxn since 

2 



\x-x r = 



x - 2 ~ X o 



\{X Q +X^ 



,WX G SSK nxn , 



(for example, see for details). 



Note that the structures of the general solution and a least squares solution of 
the matrix equation AXB — C when the equation is consistent and inconsistent, 
respectively, are exactly the same. Therefore, the structures of the solutions 
of Problems Q] and [5] are the same, too. The former and the latter facts are 
immediately seen from Theorem [T] and [21 respectively, together with Lemma [3] 
in the framework of Lemmas [1] and [H 

We conclude the paper by giving a few numerical examples. As it was 
mentioned earlier, the matrices in the examples have been taken from related 
reference work articles. It is seen that the solutions are exactly the same with 
four decimal digits as those in the works cited. All the computations have been 
performed using Matlab 7.5. 

Example 1 U 11 Example 4]- Consider the skew-symmetric solution of Prob- 
lem^ where 



A 



1 


3 


-5 


7 


-9 


2 





4 


G 


-1 





-2 


9 


6 


-8 


3 


G 


2 


27 


-13 


-5 


5 


-22 


-1 


-11 


8 


4 


-6 


-9 


-19 



B 



4 





8 


-5 


4 


-1 


5 





-2 


3 


4 


-1 





2 


5 





3 


9 


2 


-6 


-2 


7 


-8 


1 


11 



7 



c = 



171 


-537 


74 


-29 


-281 


142 


-278 


212 


-92 


-150 


196 


-523 


-59 


-111 


24 


661 


-1507 


922 


-234 


-1003 


-39 


-192 


-207 


186 


-227 


-165 


-292 


-1154 


76 


422 



j x — 



1 





4 


-1 





5 


3 


2 


7 


4 


-1 


-2 





-1 





2 


6 


1 


8 


-4 





3 


1 


4 


2 



By the formula \4-?ty in the framework Remark^ the skew-symmetric solution 
matrix is obtained as 



X = 



0.0000 2.0000 

-2.0000 0.0000 

1.0000 -2.0000 

2.0000 -1.0000 

0.0000 4.0000 



-1.0000 -2.0000 0.0000 

2.0000 1.0000 -4.0000 

0.0000 -1.0000 -0.0000 

1.0000 0.0000 -4.0000 

-0.0000 4.0000 0.0000 



Example 2 \14\ Example 
where 



1]. Consider the symmetric solution of Problem^ 



A = 



E55 
Z45 



Z54 
Pi 



B = 



C 



Z54 



Z 45 
H 5 



X 



■Z54 
h 



Z45 
Z55 

h 



Here E mn and Z mn are mxn matrices whose all entries 1 and 0, respectively, 
and P4 and H§ denote the 4x4 symmetric Pascal matrix and the 5x5 Hilbert 
matrix, respectively. Moreover, T4 and K4 are Toeplitz and Hankel matrices 
given by 



T A = 



12 3 4 

2 12 3 

3 2 12 

4 3 2 1 



id Ka = 



12 3 4 

2 3 4 

3 4 

4 



respectively. By the formula the symmetric solution matrix is obtained as 



0.8258 


-0.2692 


-0.2480 


-0.2214 


0.4129 














-0.2692 


0.6358 


-0.3430 


-0.3164 


0.3179 














-0.2480 


-0.3430 


0.6783 


-0.2952 


0.3391 














-0.2214 


-0.3164 


-0.2952 


0.7314 


0.3657 














0.4129 


0.3179 


0.3391 


0.3657 


1 





























1 





























1 





























1 





























1 



Example 3 f PA Example 3.1]. Consider the symmetric solution of Problem^ 
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where 



A 



4 


3 


-1 


3 


1 


-3 


2 


3 


-2 


3 


-4 


3 


2 


1 


4 


3 


-1 


3 


1 


-3 


2 


3 


-1 


3 


-1 


3 


2 


1 


4 


3 


-1 


3 


1 


-3 


2 


3 


-1 


3 


-1 


3 


2 


1 



-3 


4 


-3 


-3 


4 


4 


5 


-3 


5 


5 


-3 


-3 


-6 


2 


-6 


-6 


2 


2 


do 


4 


-8 


-8 


4 


4 


4 


-5 


4 


3 


-2 


-7 


-3 


2 


-3 


-3 


2 


2 


-1 


-2 


-1 


-1 


-2 


-2 



43 


-54 


73 


-54 


51 


-54 


-31 


37 


-61 


37 


-53 


37 


43 


-54 


73 


-54 


51 


-54 


-31 


37 


-61 


37 


-53 


37 


47 


-54 


73 


-54 


21 


-54 


-31 


27 


-61 


27 


-53 


27 



, X 



-1 


2 


-3 


2 


-1 


1 


00 


2 


-1 


3 


-3 


2 


-3 


4 


-3 


3 


-3 


3 


-2 


1 


-1 


2 


-3 


3 


2 


2 


2 


4 


-3 


2 


-2 


2 


-1 


3 


-3 


4 


3 


1 


1 


2 


1 


1 


1 


-2 


1 


3 


4 


-1 


1 



By the formula j4-l\ l in the framework Remark^ the symmetric solution matrix 
is obtained as 



1.7699 


1.8581 


-3.5455 


2.8924 


0.5920 


0.8523 


2.3693 


1.8581 


-0.6722 


1.8908 


-1.9156 


3.1173 


0.5698 


-1.0561 


-3.5455 


1.8908 


-0.5812 


1.3562 


-3.9861 


1.7472 


2.2490 


2.8924 


-1.9156 


1.3562 


-3.8543 


-0.6224 


0.3241 


3.6833 


0.5920 


3.1173 


-3.9861 


-0.6224 


-2.3618 


-2.2044 


3.2716 


0.8523 


0.5698 


1.7472 


0.3241 


-2.2044 


-0.0556 


2.7992 


2.3693 


-1.0561 


2.2490 


3.6833 


3.2716 


2.7992 


0.0308 
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